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ABSTRACT 



The thermal signature of a surface vehicle, "SPRUANCE” 
type destroyer, quantified by the signal to noise ratio of a 
thermal seeking missile head is parametrically investigated. 
The effects of the ship internal temperature, ship body 
paint emissivity, sky condition, sun elevation angle above 
the horizon, atmospheric profile, missile optics and flight 
altitude are examined in detail. Results show that both the 
ship body temperature and signal to noise ratio increase as 
the incident solar energy and the ship body paint emissivity 
increase and that the signal to noise ratio appears a peak 
for sun elevation angle in the range of 40° to 60°. Moreover 
the signal to noise ratio increases as the missile flight 
altitude decreases and keeping the other parameters constant 
higher values are found for the Midlatitude Summer atmos- 
pheric profile than for the Tropical atmosphere. 
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I. INTRODUCTION 



The basic function of a thermal-seeking-head missile is 
to track on a surface vehicle using the vehicle's radiant 
intensity. This operation is based on the fact that a moving 
vehicle shows a generally higher surface temperature that 
the environment. Reasons for this are the existence of 
operating engines, resulting in both increased engine room 
temperature and ''hot” combustion products exiting into the 
atmosphere from the stack, the operation of a number of 
devices which result in heat generation, the requirement for 
certain compartment temperatures comfortable to the crew, 
and the incident solar energy on the painted metallic 
vehicle surfaces . 

Exact information about the heat seeking missile charac- 
teristics is not available because of their classification. 
Hence our study will be a parametric one in order to examine 
the relative influence of a number of parameters on the 
vehicle thermal signature which is quantified by the ’signal 
to noise ratio of the thermal detecor head. The range of the 
parameter values will be chosen using both experience and 
unclassified information. 

Exceptional cases such as operation of the vehicle close 
to the coast or other vehicles, which could have similar or 
higher temperatures, will not be considered. At this point 
we must realize that the thermal interaction of the vehicle 
with the environment is a very complicated phenomenon, and 
its exact determination for the different parts of the ship 
body is outside of the scope of the present work. Instead we 
will consider a simplified case of uniform interaction for 
all the ship parts with the environment (convection, radia- 
tion, solar absorption energy) and will limit the source for 
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the thermal head detector enviromental noise only to that of 
the sea water. 

The main procedure steps of the present study will be : 

To create a simplified model for the calculation of the 
vehicle body temperature and to show what parameters affect 
it and in what direction. 

To determine the "hottest” vehicle area and combining 
with the operating characteristics of the missile to calcu- 
late the dimensions of the tracked area of the vehicle. 

To calculate the total radiant intensity of the tracked 
area of the vehicle in a certain wavelength range, where the 
thermal detector head is operating. 

To calculate the atmospheric transmittance, using 
LOWTRAN 6, for each given combination of the parameters of 
our interest in the given detector wavelength range and 
consequently the radiant heat flux detected by the thermal 
head . 

To create an appropriate detector model that could show 
optimum operation for the conditions of our problem, and 
finally to calculate the creating signal- to-noise ratio 
which quantifies the vehicle thermal signature and describes 
the detector operational quality. 

To discuss the relative influence of all the pertinent 
parameters and to propose certain considerations about both 
the vehicle and the missile. 

Without expecting significant quantitative importance 
for our results, due to the simplifications that have been 
made, their qualitative importance must be more remarkable. 

Parameters that will play significant roles in our study 
will be : 

vehicle internal temperature, 

vehicle body paint emissivity, 

sky condition, 

sun elevation angle above the horizon, 
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created by the 



atmospheric profile, 
missile flight altitude, 
instantaneous Field of View (I.F.V.) 
detector's optic system, 
detector material, 

The present work consists of five Chapters. In the first 
Chapter we state the assumptions of the problem, we calcu- 
late the vehicle body temperature, and we determine its 
"hottest" area which is tracked by the detector. In the 
second Chapter we use the LOWTRAN 6 program in order to 
determine the atmospheric transmittance and in the third 
Chapter we calculate the total vehicle and background 
radiant intensity incident on the tracker aperture. In the 
fourth Chapter we calculate the detector signal- to-noise 
ratio for all the possible combinations of the parameters , 
and finally in the fifth Chapter we discuss the relative 
importance of each parameter and we suggest ways that could 
improve both the vehicle thermal signature and the detector 
operation. 
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II. SHIP MODELING 



In order to create an exact ship model, it is necessary 
to acquire a significant amount of empirical data. For the 
present study, and for several reasons, the available data 
are limited, and are based only on the unclassified litera- 
ture given in the first six references, with relatively 
reasonable assumptions. 

A. ASSUMPTIONS 

All compartments of the ship have a uniform temperature 
of 20°C=293.15K except the main engine rooms 1 and 2, which 
have a temperature of 30 0 C= 303 . 15K . 

The stack exit temperature is uniform and we study the 
heat radiation from the exit plane neglecting the plume 
radiation. Reasons for doing this are the non-uniformity of 
the plume geometry and temperature distribution, and the 
fact that it is cooled significantly in a relatively short 
distance from the exit plane [Ref. 3]. We choose as stack 
exit temperature 306 0 F= 152 . 2 0 C=425 . 4K 

In order to account for the heat convection from the 
ship to the air we will use the relative wind-speed of 30 
Knots. For this case we will have forced convection with a 
heat transfer coefficient in the range of 5 to 10 (BTU/H 
FT 2 R) or equivalently 28.4 to 56.8 (W/M 2 K). We will not 
account for the heat leakage from the ship body to the sea 
water or for the heat loss due to the vaporization effect on 
the wet surfaces above sea level. For this reason and in 
order to be closer to reality we choose for our heat convec- 
tion losses calculation the value of ID (BTU/H FT 2 R) or 56.8 
(W/M 2 K) as representative for the convection heat transfer 
coefficient . 
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The present study will be in Mediterranean atmospheric 
conditions and more specifically for the area of the Aegean 
Sea. We will use two atmospheric profiles i.e. Midlatitude 
Summer for months October through April and Tropical for 
months May through September. The corresponding ambient 
temperatures are 21°C=294.2K and 26 . 6 0 C=299 . 7K respectively 
[Ref. 6]. 

B. GEOMETRY AND MATERIAL OF THE SHIP 

In Figures 2.1 and 2.2 you will find the ship side and 
top view respectively with the appropriate dimensions. 

The material of the ship lower body construction is 
constructional carbon steel with thickness 13.5mm and 
thermal conductivity 54 (W/M K) [Ref. 9]. 

The material of ship body superstructure is duralaminum 
with thickness 3mm and thermal conductivity 164 (W/M K). 

These small values of the thickness and the high values 
of the thermal conductivities for both materials result in 
negligible temperature gradients inside the materials. Hence 
it is reasonable to assume a uniform temperature throughout 
the thickness of both materials. 

Both structures are painted with gray paint i.e. the 
emissivity and absorptivity will be spectrally independent 
and equal between them. A typical value for oil paint emis- 
sivity - absorptivity is 0.94 [Ref. 10:p.07] and for TiQ 
gray paint is 0.87 [Ref. 4:p.278] both for temperatures 
300K. From these values we observe that the available data 
give us paint emissivity - absorptivity values in the range 
of 0.90. Based on the facts that these values play a signif- 
icant role in the determination of the ship body temperature 
and that intensive research resulted in the production of 
new paints with half the values of the already existing ones 
[Ref. 7] we will use two (2) different values i.e. 0.80, 
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0.55. Therefore we cover the useful range of the existing 
and the newly developed paints. At this point we must note 
that we don't have any available information .about selec- 
tively absorbing paints that could exist and give the 
minimum possible thermal signature and our assumption about 
the equality between emissivity and absorptivity is inevi- 
table as a first approximation. Moreover the solution of the 
problem for two (2) different emissivities constitutes a 
good research point for the relative influence of this 
factor . 

The emissivity of the interior surface of the ship body 
will be approximated as that of stainless steel which is 
0.44 [Ref. 4:p.278j. 



C. SOLAR HEAT FLUX IN THE SEA LEVEL SURFACE 

In Figure 2.3 you will find the illuminance levels on 
the surface of the earth due to the sun for three different 
sky conditions i.e. Unobscured Sun, Sun with Light Clouds 
and Sun with Heavy Storm Clouds as a function of the eleva- 
tion angle above the horizon. These three sky conditions and 
the elevation angle above the horizon will be two more 
parameters in the solution of our problem. 

Since the given values in Figure 2.3 are in photometric 
units (Lux) and we want to convert into useful radiometric 
units (W/M 2 ) we have to determine a conversion factor. 

We know that the maximum incident solar flux on the 
earth surface is approximately 900 (W/M 2 ) [Ref. 10:p.ll]. 
This value coresponds to 1.15* 10 s (Lux) and therefore the 
conversion factor will be: 

CF = 900 / 1 . 15 • 10 5 = 5.826* 10' 3 (W/M 2 Lux ) 
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Figure 2.1 Ship Side View 
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Figure 2.2 Ship Top View 
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Using this conversion factor and the values of Figure 
2.3 we obtain the normal component of the solar heat flux in 
(W/M 2 ) on a surface at sea level. Table 1 demonstrates the 
above values . 

In order to calculate the values of the solar heat flux 
falling obliquely on a surface we must introduce the inci- 
dent angle 0 and multiply the values of Table 1 by cos0. 

For a vertical face the angle 0 is identical with the 
elevation angle above the horizon. For a horizontal face 
i.e. the deck in our case, the angle 0 is the complement of 
elevation angle above the horizon. 

D. CALCULATION OF SHIP BODY TEMPERATURE 




Figure 2.4 Control Volume for Determination of Ship 

Body Temperature 

In Figure 2.4 you will find the appropriate control 
volume of the ship body and the pertinent energy flows that 
are related with it. For simplicity the control volume has 
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unit area exposed to both the environment and the ship 
interior. More specifically: 

6* : paint emissivity 

a„ : paint absorptivity (£<►= a<, ) 

: emissivity of the interior surface. ( £«= 0.44 ) 

T 0 : ambient temperature 
T,», : ship compartment temperature 
E w : normal component given by Table 2 
Q 0 : net heat exchange of control volume with the 
environment 

Q, n : net heat exchange of control volume with the 
interior of the ship 

h : convection heat transfer coefficient with the air 
The absorbed solar heat flux is E^a©. The net heat 
exchange with the environment is due to heat radiation and 
heat convection i.e. 



The net heat exchange with the interior is due only to heat 
radiation i.e. 



Energy balance of the control volume under thermal equilib- 
rium gives : 



Q e =£c cr( T 4 -T 4 )+h( T-T. ) 



( 2 . 1 ) 



Qm=^( T 4 -T 4 n ) 



( 2 . 2 ) 



E sn «a *=£.cr( T u - To ) + h(T-T.)+Qr( T 4 - T,J, ) 



(2.3) 



or 



F ( T ) = (e.+ e, n ) crT 4 + hT - <7 (e, T. 4 + ) - hT e - E sn « a 0 = 0 



(2.4) 
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Equation 2.4 must be solved numerically. We will use the 
Steffensen Recursive Method [Ref. 11]. The idea of this 
method is that starting with an x close enough to the exact 
solution S of the equation F(x)=0 then the sequence: 

x M fx K - [f 2 (x K )/ f (x*+f (x K ) )-f (x K )] (2.5) 

will converge very fast to the exact solution S. 

In our case the equation is F(T)=0 and we are looking 
for the value of T that satisfies it. 

In order to obtain an initial value so that our sequence 
will converge we use the value: 

Tin" [ ( E sn * a, + *.°V + 4-PTm ) / (^ + ^) <T] u « (2.6) 

The above initial value of T does not account for the 
heat convection and for our purpose gives acceptable initial 
values. Finally we stop the recursive formula when |T U1 ,-T K | 
is less than or equal to 10* 6 . 

In Appendix A, Tables 5 through 18 demonstrates the ship 
body temperature for both horizontal and vertical faces, for 
all the combinations of two (2) atmospheric profiles, three 
(3) different sky conditions, two (2) internal (T, n ) tempera- 
tures and two (2) paint emissivities . 

E. DETERMINATION OF THE SHIP "HOTTEST" AREA 

From the calculations in section D, we observe that the 
higher ship body temperatures happen to be in the region of 
the two engine rooms (higher internal temperature) . 

From the ship geometry^ it is obvious that the dimensions 
of the engine room No 1 are greater than those of engine 
room No 2. Therefore the total radiant intensity (W/Sr) of 
engine room No 1 will be greater since for both regions we 
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have the same normal components of radiance (W/cm 2 Sr). In 
addition to that to each engine room there corresponds the 
same stack geometry also going a 'hot' region too. 

Under this simplified assumption we conclude that for a 
heat seeking missile head the area of engine room No 1 from 
the stack exit all the way down to sea level is the most 
significant region. 

At this point it is necessary to couple the ship with 
some characteristics of the missile detector operation. 

Without getting involved in details of the detector 
operation, which will be studied in later sections, we will 
determine the detector Instantaneous Field of View (I.F.V.), 
so that we will know exactly the dimensions of the 'hottest' 
area which will be tracked by the detector and therefore we 
will be able to calculate the amount of radiant energy that 
will be transmitted to the missile. 

F. DETERMINATION OF THE DETECTOR INSTANTANEOUS FIELD OF 

VIEW 

The following range of data is based on our experience 
and research, aiming to calculate a reasonable and useful 
I.F.V. In Figure 2.5 you will find the corresponding 
geometry . 

The range of the missile flight altitude will be from 
50m to 100m above sea level. Since we want to reduce as much 
as possible the background noise, the detector must have an 
inclination with respect to the horizontal. We choose the 
range for this inclination (angle ) from 2° degrees to 5° 
degrees . 

From general cosideration for the detectors [Ref. 4] we 
know that I.F.V. is a function of the 'hottest' ship area 
that the detector has to track. 
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For the creation of I.F.V. it is necessary that the 
detector accepts radiation from inside an area which is 
related to an angle 25 generated by the detector mirror 
system as in Figure 2.5 . From the same figure we conclude 
the following: 

8 = 90°-<f -6 
BC = AB + tan0 
AC = (AB 2 + BC 2 ) 1,1 
CF = AC*sin6 

CG = 2CF = 2AC»sinS = SIDE 
AF = CG/2tan6 = PATH 
Where : 

CG represents the side of the I.F.V. 

AF represents the distance or path length between detector 
and the I.F.V. 

AB = H : flight altitude 

For simplicity we will consider that the shape of I.F.V. 
is square with side equal to CG of Figure 2.5 

In Appendix B Tables 19 through 29 demonstrate the 
values of PATH and SIDE for all the possible combinations of 
S and H. Program SHIP 1, in Appendix E, has been used for 
the calculation of the above values. 
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Determination of Instantaneous Field of View 



G. DETERMINATION OF THE DIMENSIONS OF THE TRACKED SHIP AREA 



In the previous section we have shown that the detector 
will track on the area of engine room No 1. 

In Figure 2.6 you will find a three-dimensional sketch 
for the engine room No 1 area. Before we start the calcula- 
tions of the different surfaces we will explain why the side 
of I.F.V. is greater than the width of the engine room No 1 
which is Wi=l, 620cm. We consider as useful that the I.F.V. 
will include the whole area between sea level and the stack 
exit plane. The total height of this dimension is: 

^ + H 3 + H<,+ H & =720+900+360+180=2, 160 cm 
From the Tables 13 to 24, we choose for each flight 
altitude and for constant 5=0.56° the minimum SIDE that is 
greater than or equal to 2.160cm 

Table 2 demonstrates the results of this choice 
Based on Figure 2.6 we will now determine the character- 
istics of each individual surface. 

A x = Wi • Hi 
Vertical face 
T, n = 303. 15K 

4r= (SIDE - Wi )*Hi 
Vertical face 
T, n = 293. 15K 

A a. = W^ ^ 

Horizontal face 
T,„= 303. 15K 

• 

A n = (SIDE - W x )*Hi 
Horizontal face 
T,*= 293. 15K 
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Figure 2.6 Three-Dimensional Sketch of the "Hottest" Area 



TABLE 2 

Final Choice for the Dimensions of I.F.V. 



HEIGHT 

FLIGHT 


SIDE 

(m) 


PATH 

(m) 


DELTA 


PHI 


50 


(m) 


21.88 


1119 . 34 


0.56° 


2.0° 


60 


(m) 


21.97 


1123.89 


0.56°. 


2.5° 


70 


(m) 


22.04 


1127.25 


0.56° 


3.0° 


80 


(m) 


22.09 


1129.84 


0.56° 


3.5° 


90 


(m) 


22.13 


1131.96 


0.56° 


4.0° 


100 


(m) 


22.16 


1133.72 


0.56° 


4.5° 



Aj = SIDE • H 5 
Vertical face 
T in = 293. 15K 

A s = Wn» H* 

Vertical face 
T, n = 303. 15K 

At this point we must notice that surface A s corresponds 
to the stack and, since we don't have information about the 
insulation of this structure, we approximate its temperature 
with that of the surface with T = 303.15°K 
A* = 4 . 43 • 10 “* cm 2 
Horizontal face 
T, n = 425. 4K 

A*= SIDE-H*- 0/4)^ 

Horizontal face 
T 1V> = 29 3. 15K 
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For the order of magnitude of PATH i.e. the distance 
between missile and ship we will approximate the thermal 
emission of each individual surface as point source emis- 
sion. These "point sources" will be: 

(A x + An) , (Ai + Aji) , A 4 , At, , A s , A e , A f . 

with corresponding mean altitudes from sea level: 

H*i= H Hj. / 2 

h**= ^ 

H« 4 = H t + H 5 / 2 
Hw,- Ht + H 4 + H„ / 2 
H» % = H A + H 5 + Ht, + H s /2 
Hmq = H t + H 4 + + Ht, 

H m ,= H a + Hj, 

The data for the above calculations are: 

Wj. = 1,620cm 
Hi = 720cm 
H x = 210cm 
Hj = 900cm 
H^ = 1,260cm 
Wj, = 900cm 
H^ = 360cm 
H 5 = 180cm 
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III. ATMOSPHERIC TRANSMITTANCE USING LOWTRAN 6 



A. DESCRIPTION OF LOWTRAN 6 

LOWTRAN 6 is the latest issue of LOWTRAN code that has 
been designed to calculate the transmittance of the atmos- 
phere over spectral bands from 0.25pm to 28. 5pm along paths 
of various inclinations and, of course, through different 
atmospheric profiles and parameters. 

In order to use LOWTRAN 6 program we have to fill four 
(4) cards (groups) of input data. 

Cards 1 and 2 require the significant meteorological 
parameters for the atmospheric profile. Card 3 requires the 
geometry of the problem i.e. initial and final altitude and 
the path length. Card 4 is related to the wavelength range 
and limits . 

B. DATA USED 

For card 1 we used data for the Mediterranean Sea and 
more specifically for the Aegean Sea area and the region 
from Crete to North Africa. 

For this area we have two different atmospheric 
profiles, namely Midlatitude Summer for months October 
through April, and Tropical for months May through September 
[Ref. 6]. 

For card 2 we used again data from [Ref. 6] for each 
month individually and for card 3 we used all the possible 
combinations between final altitude and PATH from Table 2 
and for the initial altitude the mean altitudes of each 
"point source" of the "hottest" ship area as given in 
Section G of Chapter II. 
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At this point we have to choose one more parameter that 
will be needed for our problem : the wavelength range of 
the detector operation. We decide that this range will be, 
in first approximation, in the 5um window (i.e. A\=4.20pm - 
>^=5.18^111) [Ref. 4 : ch 5, p.89]. 

Lowtran 6 has the advantange of giving us, besides the 
average transmittance between two limits, the spectral 
transmittance as well. From this spectral distribution we 
can reject certain values of very low spectral transmittance 
and therefore the wavelength range can be narrower. This 
peculiarity will have an influence on card 4. 

C. RESULTS 

Running LOWTRAN 6 for all the flight altitudes of the 
missile and for both Midlatitude and Tropical atmospheric 
profiles, we observed that, for wavelengths less than 
4.435pm and greater than 5.0pm, the spectral transmittance 
generaly shows values less than about 0.10. Deciding to 
reject these low values we choose as the final spectral 
range of the detector operation to be from 4.435iim to 
5 . Oum 

In Figure 3.1 you will find the spectral transmittance 
in the above range for both atmospheric profiles. 

Using the final wavelength range we calculated the 
average transmittance for all the possible combinations 
between atmospheric profile, missile flight altitude and 
individual "point sources" of the ship's tracked area. 
Tables 23 to 46 in Appendix C demonstrate both the data and 
total average transmittance for each month. 

From the data of the above tables, it is obvious that 
the influences of both missile flight altitude (from 50m to 
100m) and the relative differences between the individual 
"point source" altitudes (4m to 16m) are very weak on the 
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average transmittance. Averaging in both directions in this 
two dimensional table results in a unique representative 
value for each month. Checking these average values we 
observe again a relatively small difference between months 
belonging in the same atmospheric profile. Since a unique 
value of average transmittance of each atmospheric profile 
would be very convenient for the calculations to follow, we 
average over several months yielding the following final 
values : 

Midlatitude summer : 0.4711 
Tropical: 0.4259 
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WAVELENGTH IN MICRONS 



IV. CALCULATION OF SHIP AND BACKGROUND RADIANCE 



A. NORMAL COMPONENT OF SHIP RADIANCE DUE TO THERMAL 
RADIATION 

. The purpose of this calculation is to find, by inte- 
grating the Plank's radiation law, the radiance in the given 
wavelength range. 

From Plank's radiation law we know that: 



c* =2TTc 2 h=3.7Al5‘10 u (Wcm 2 ) 

Cx =hc/k= 1 . 4388 (cm*K) 

The above formula is in the form of: 

M (T) = C*x <Tn (£* -1) 

Integrating this between the band limits , we will 

find the in-band energy radiance [Ref. 4:p. 19]. 

Therefore : 



M^T) = [2frc 2 h/ X s ] * [1/ (e K>>T -1)] (W/cm 2 cm) 



or 



M (T) =C^ 5 /(e^ - 1) 



where : 



Where : 






and x=Ci /aT 



CO 



21i=2Tm* ,i *m''* [ (m»x ) 3 + 3(m*x ) 2 + 6m*x+6] - 
e***m 4 [(m*x ) 3 +3(m # x ) 2 + 6m*x+6] 



In the above relations^ the temperature T 
Tables 5, 18 in Appendix A. 



is taken from 



B. NORMAL COMPONENT OF SHIP RADIANCE DUE TO REFLECTED SOLAR 



ENERGY 

We have already calculated the normal component of solar 
irradiance E incident on each ship surface. In the energy 
balance calculation we considered that the solar energy 
absorbed by the surface is l^a,. Obviously- the amount 
(l-a 0 )»Ej n is reflected. 

Defining as D(0-> ) the percentage of the solar constant 
associated with wavelengths shorter than A we take 
[Ref. 4 : p . 17]. 

D(0-4. 435) =99. 3075454% 

D(0-5.000)=99. 511500% 

Therefore 

AD=D (0- 5. 000 )-D( 0-4. 435) =0.2039546% 

In conclusion we have that the normal reflected solar 
energy per m 2 from >» = 4. 435pm to \t=5.00pm is: 

(1-a* ).AD*E sn =(l-a t> ) • 2 . 04* 10* 3 *E** (W/m 2 ) 

In order to determine the normal component of radiance 
due to reflection we must divide the above value by tr and in 
order to transform into* cm 2 we must divide by 10“ . Finally 
the normal component of radiance due to reflection in 
(W/cm 2, Sr) from Xi = 4. 435pm to >i=5.00pm will be: 

Lra>= ( 1 - a 0 ) • 2 . 04* 10' 3 ‘E^/ir • 10 4 = ( l-a 0 ) »4 . 08 • 10' 8 »E 30 (W/cm 2 »Sr) 

C. NORMAL COMPONENT OF STACK EXIT PLANE RADIANCE 

So far we know that the stack exit plane area is 
4. 43 *10 4 cm 2 and that the corresponding uniform temperature 
is 425.4K. 

For simplicity we consider that the greater part of the 
exaust gas is CO* . Since COj. is a selective emitter we need 
to carry out the approximation shown in Figure 4.1 which 
partially accounts for atmospheric absorption and approxi- 
mates the peak with a constant emissivity over a limited 
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Figure 4.1 The 4 . 4um Emission Band of COj[Ref. 12] 

range. Taking the band to be from %=4.35um to ^i=4.55um 
and the emissivity 0.5 we find for the radiance: 

L =eo/M>dA = ( €»ctT 4 /tt ) / M»dW/ M x dx= (e.crT'*/^)* (0.05) 

-iu}5 4i5 ■'» 

= (0.5)* (5.67* 10* 11 )• (425.4) U# (0.05)/Tr 
= 1.478* 10' 3 (W/cm 2 *Sr) 

Since our band of interest overlaps with half of the 
4.4pm emission band of CO* the above value must be divided 
by 2. Therefore: 

L=7. 39*10^ (W/cm 2 *Sr) 

D. TOTAL SHIP RADIANT INTENSITY IN THE DIRECTION FROM SHIP 
TO MISSILE 

As mentioned in Section G of Chapter II, we will consider 
the surfaces (A a +A U ) , (Ax+Aj,), A s , Ai, , A 5 , A^ , A,, as "point 

sources". For each source we know the corresponding radiance 
due to thermal radiation and reflection i.e. 

Ln =L&>, +Ltuu v 
L-tu =L/un + L^ 4 ^„ 

Ltx =L^+L«.&x x 
Lt« =Ln A x + Lt U i.« 

las = Ltos + LR£,\ 5 
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Lta, =L,6Xl, + L>Rfc>o< 

Lts = L^ s + L?.&xe, 

L-t? = Lt,yx + Le > *\> 

For L-rt we have the unique value given in section C of 
chapter III 

In order to obtain the normal radiant intensity (W/Sr) 
of each surface we must multiply the corresponding values of 
each area and radiance. Doing this we take: 

I iu =Ln •A l + L™ * A V1 
I in - Ln * A* + Aii_ 

I»h = Lt»*Aj 
Ii*h = Lri»*A i , 

Is*» = Lts*A B 

I«H = Lt 4*A 6 
I**» =Ln*A> 

The important magnitude for us is the radiant intensity 
in the direction of source - missile. 




Figure 4.2 Relative Position Between Surface- Missile 

In Figure 4.2 you will find a sketch representing the 
relative position between the surface and the missile. Based 
on that we conclude this the radiant intensity will be: 
I=In*cosx for a vertical face 
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I = If**sinx for a horizontal face 
where x=sin 1 (H-Hia/PATH) 

The resulting relations will be: 

Ij. =Iirt • cosx 
I a, = I»i* sinx 
It, =IjfV* cosx 
I<, =1^* cosx 
I 5 =I 5I »* cosx 
1^ =ItH* sinx 
I* = !,,»• sinx 



E. BACKGROUND RADIANT INTENSITY 

In the ideal case in which the detector tracks perfectly 
on the ships "hottest" area and the square I.F.V. will fit 
perfectly on that the detector will still see some areas 
outside of the ship. This is due to the fact that the upper 
construction of the ship has not a perfect rectangular 
shape. Therefore the detector will receive signal from an 
area of the sea that radiates behind the ship. Since this 
signal constitutes noise for the detector, it has important 
meaning . 

In Figure 4.3 you will find the geometry from which we 
will determine the background radiant intensity. We must 
notice that the sketch in Figure 4.3 has been magnified for 
convenience in the ship region. In other words angle (0+S) 
is very close to 90° and the side of I.F.V. AE is ASE is 
practically equal to the ship dimension ABC. 

AB represents the height of surface A^ [Figure 2.6] from 
the sea level or AB=H i +H X 

BC then represents the path through which' the detector 
will see the sea area FG past the stack. 
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Figure 4.3 Geometry for Determination of Background 

Radiant Intensity 



The width of this area is (SIDE-W A ) which represents the 
free horizontal dimension over the surface A*. 

From the geometry of Figure 4.3 we have: 

DS=AS- SIDE 
0X =tan' 1 (DS/PATH) 

0 » = 6 - 0 !. 

MA=PATH/cos6 
0= cos' 1 (H/MA) 

PF=H»tan(0<-6 + 0i ) 

PG=H* tan(0+ 26 ) 

FG=PG-PF 

0j =0 + 5 + 0i +0x / 2 

MK = H» tan0 5 

AREA=FG* [SIDE-W 6 ] • cos0& 

Table 3 demonstrates the result of the above calcula- 
tions for six different flight altitudes. 



TABLE 3 



Geometrical Parameters for Observed Sea Area 



H (m) 


FG (m) 


0 (deg) 


MK (m) 


AREA (it?) 


50 


333.67 


88.415 


1806.98 


144 . 90 


60 


230.54 


87.915 


1648.07 


131.68 


| 70 


180.09 


87.410 


1547.47 


127.76 


f 1 

80 


144.50 


86.91 


1481 . 9*-. 


122.29 

• 


, 90 


116.07 


86.4 


1430.50 


114.42 


: 100 


9 a : 4 : 


85.9 


1395 . 07 


111.60 



For the sea water we know that the average temperature 
in the geographical area of our interest is 13°C=286.15K for 
Midlatitude Summer and 18 ° C= 29 1 . 15Kf or Tropical atmosphere 
[Ref. 6]. 
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For the sea water emissivity we chose the value of 0.96 
[Ref. 10 : p . 07] . 

For the calculation of the normal component of the sea 
radiance (W/cm 2 *Sr) in the range !>u=4. 435pm to ^=5. 00pm we 
will use again the formulas from section A of chapter IV. 
This value must be multiplied by the water emissivity 
=0.96. For convenience we consider the observed area of sea 
as a point source. 

The normal component of the radiant intensity of this 
point source it will be: 

Iw* = L w jAREA 

And in the direction leading to the missile will be: 

I* = I -H » cos© 
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V. DETECTOR MODELING AND CALCULATION OF S/N RATIO 



A. DETECTOR SYNTHESIS 

We already studied the geometry of the detection problem 
which shows us what the detector is required to do, and of 
course from this side we are interested in the most real- 
istic performance of the detector. Since we do not have 
available experimental data that could show us the perform- 
ance of the detector element we have to choose this 
according to the literature. 

From Figure 5.1 [Ref. 10] the possible choices for the 
detector include photovoltaic or photocontactive Indium 
Antimonide and Lead Selenide. The detailed data sheets for 
these detectors [Ref. 12 : pp . 3 65 - 3 68] show the relative 
advantages and disadvantages of each. When they are compared 
on the basis of D' v , photovoltaic Indium Antimonide is 
clearly the best detector. A possible difficulty with this 
conclusion arises when the responsivities are compared and 
it is noted that value of Lead Selenide is 10 times higher 
than it is for the other two types of detector. In an ideal 
operating environment, that is, in which there are no extra- 
neous sources of noise, the designer is interested in the 
signal - to - noise ratio of the detector, as described by 
D" and not the signal level, as described by responsivity . 
In advance we know that during the experimental trials 
photovoltaic Indium Antimonide had the better performance. 
Therefore we are led to the decision to use one element of 
photovoltaic Indium Antimonide, and we can see the detailed 
data for this material in Figure 5.2 [Ref. 12:p.368]. 

This one element will be placed • in the middle of an 
assumed two mirror system as in Figure 5.3, with 
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Figure 5.1 The D* of the Detectors in Each Atmospheric 
Window at the Most Common Operation Temperature [Rex. 10] 
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Indium Antimonide (InSb) 

Type: Photovoltaic. 

Element From 0. 1 x 0. 1 to 10 x 10 mm. 
sizes Rectangular elements from 

available: 0.1 mm wide. Circular elements 

with diameters from 0. 1 to 10 mm. 
Less than 1 y. sec. 
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Figure 5.2 Data Sheet for Photovoltaic Indium 
Antimonide Detector [Ref. 12] 
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reflectivity 0.9 each mirror and total of two ^ 
=0 . 9 • 0 . 9=0 . 8 1 . The detector will be considered to be used 
with a two-mirror system designed so that the detector 
element will subtend area equal to the I.F.V. at a distance 
of PATH. Thus we have to determine the size of the sensitive 
element (detector). 




Figure 5.3 System of Two Mirrors in the Receiver 

We know that the area of the detector A is given by: 

A d = w»f 2 

where 

w=Instantaneous Field of View in (Sr) 

w=( SIDE •SIDE)/ (PATH) 2 =(2200 *2200)/ (PATH) 2 (Sr) 
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f=equivelant focal length 

and for PATH we have to use the values from Table 2 which is 
function of the flight altitude. 

For the our purpose of work a very good and accurate 
approach to calculate the equivalent focal length of the 
system of two mirrors is the following. 

First we have to choose the diameter of the entrance 
aperture D 0 =12cm. This dimension is good enough to have an 
aerodynamic shape in the nose of the missile. Thus from the 
Figure 5.3 we have: 

(AB) 2 = (AO) 2 + (OB) 2 =72 and AB=8.4cm 

(AO) 2 = (AC) 2 + (CO) 2 or 

(CO) 2 = (AO) 2 - (AC) 2 =(A0) 2 - (AB/2) 2 = 

= (CO) 2 =6 2 -(8 .49/2) 2 = 17 .98cm 2 
and C0=4.24cm 

We can define the distance CO as the equivalent focal length 

f =4 . 24cm 

And 

A<i = wf 2 = (2200) 2 ‘(4.24) 2 / (PATH ) 2 
Substituting the values of PATH for each flight altitude 

gives A and the square root of this gives us the side of 

r 'f 

the detector element .'Tab le 4 demonstrates the dimensions of 
the detector element shows that with reasonable detector 
size the I.F.V. at about 1125 meters slant path range will 
include the hot region of the ship regardless of the flight 
altitude . 



B. SIGNAL VOLTAGE 

The general formula that gives the 
V 5 = [A»* v R/ ( PATH ) 2 ] • [ J I^TqlV) dA] 

Where ** 

A 0 :Area of entrance aperture 
A 0 = TT D 2 / 4= 113 . 1cm 2 
r 0 :Optics reflectance 



signal voltage is: 
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TABLE 4 

Dimensions of the Detector Element 



FLIGHT 
ALTITUDE 
(m ) 


PATH 

(m) 


DIMENSION OF DETECTOR 
ELEMENT 

(miM) 


50 


1,120 


0.83 x 0.83 


60 


1,124 


0.83 x 0.83 


70 


1,127 


0.83 x 0.83 


80 


1,130 


0,83 x 0,83 


90 


1,132 


0.83 x 0.83 


100 


1,134 


0.83 x 0.83 



R : responsivity=2* 10 4 (V/W) 

I\ : Spectral radiant intensity of the source 
I* : Spectral atmospheric transmittance 

The above relation can be further simplified using 
instead of the spectral transmittance an average one as we 
determined in Section C of Chapter III i.e. 0.4711 for 
Midlatitude summer and 0.4259 for Tropical atmosphere. 

Taking outside of the integral then the quantity 

J I*dX becomes simply the inband radiant intensity of the 
source which has been calculated in section D of chapter 
III. Hence the signal voltage formula becomes: 

Vs=A,*r.*R*VW (PATH ) 2 

In order to determine the signal voltage due to all the 
available point sources inside the I.F.V. we use the super- 
position principle: 

Vs= [A/ r.»R •"£<*/ (PATH) 2 ]* [I k+ I 1+ I a + I li+ i s+ i t +1,] 
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C. NOISE VOLTAGE 



The noise voltage is due to two main factors , namely the 
background noise that is created from the background radiant 
intensity and the noise from the optical system. 

In Section E of Chapter III we determined the radiant 
intensity of the sea region that the detector observes and 
we named I w (W/Sr). The noise voltage of this reason will 
be : 

V N8 =A.*r.‘R*VW(MK) 2 

where MK is determined in Figure 4.3 and indicates 
the average path between the observed sea surface 
and the detector. 

The noise signal due to the optical system is given by: 
[Ref. 12]. 

V* 0 =R* (A/Af)^ /D* 

where 

D* : Detectivity with average value 9’IC) 11 [cm(Hz ) 1 / 2 /W] in 
the wavelength range ^»=4. 435pm to *^5. 00pm 

A,j : Area of the detector as determined in previous section of 
the present chapter. 

R : Respons ivity 

Af .‘equivalent noise bandwidth 
The equivalent noise bandwidth is given by 
Af=l/2Ia 

where 

Ididwell time (time required for the image of the target to 
pass across the detector) . 
and Id=w*c/£L 
where 

c:number of detector elements (1 for our case) 

Q.:Time rate of search, 
and 

Q=Q/t P 

where 
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Q:size of search field in (Sr) 

Cp: frame time (time required to scan the entire search 
field). For our case we chose 10msec. 

Therefore 

Af=l/2lj=Cl/2 , w , c = Q/2*w*c , 'Cp 



DETECTOR 




Figure 5.4 Total Detector Search Field 

In Figure 5.4 you will find the geometrical representa- 
tion of the total search field which will have total length 
of 1000m and is scanned by the detector in =10msec. 

The total solid angle that the search field subtends at 
the detector is: 

Q= dQ 

dQ=dx»SIDE/« 2 
t 2 = (PATH) 2 -x 2 
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dQ= [ SIDE / ( PATH ) 2 - x 2 ] • dx 
Q= ( SIDE ) dx/ (PATH) 2 -x 2 

Since (PATH) 2 >500 2 always, from the tables we get: 

Q= (SIDE) /2 (PATH) [ln( (PATH+x) / (PATH-x) )] 

= (SIDE) /2 (PATH) [ln( (PATH+ 500 )/ (PATH- 500 ) ) 

- ln( (PATH-500)/ (PATH+500) )] 

= ( (SIDE) /2 (PATH) )• In [(PATH+500)/ (PATH- 500)] 2 

Q= ( (SIDE)/ (PATH) )»ln[ (PATH+500)/ (PATH- 500)] 

Therefore 

Af= (1/2 (SIDE)* c*Cp) «ln[ (PATH+ 500 )/ (PATH- 500 ) ] 

The total noise voltage is 
Vn=V N% +V Ho 

And the signal to noise ratio is: 

Vs/Vn=Vs/(V N&+ V* 6 ) 

In figures 5.5 to 5.22 we can see the summary of the 
results of the above calculation for all the possible combi- 
nations of the relative parameters. 

In Appendix D we can see the influence of sky condition 
to S/N ratio for the two extremes flight altitude i.e. 50 
(m) - 100 (m) and for all the combinations between atmos- 

pheric profile and paint emissivity. 
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Figure 5.7 Influence of Missile Flight Altitude for Sun with 
Heavy Storm Clouds, £*=0.80, Midlatitude Summer 
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Figure 5.10 Influence of Missile Flight Altitude for Sun with 
Heavy Storm Clouds, £*,= 0.55, Midlatitude Summer 




62 




i 



Figure 5.12 Influence of Missile Flight Altitude for Sun with 
Light Clouds, £<>=0.80, Tropical Atmosphere 
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MISSIIF FLIGHT ALTITUDE IN METERS 
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Figure 5.16 Influence of Missile Flight Altitude for Sun with 
Heavy Storm Clouds, f»=0.55, Tropical Atmosphere 
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Figure 5.21 Influence of Emissivity and Atmospheric Profile 
for Flight Altitude 100m and Sun with Light Clouds 
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VI. RESULTS AND DISCUSSION 



A. SHIP BODY TEMPERATURE 

■Tables 5 to 18 of Appendix A show the effects of the 
ship internal temperature, atmospheric profile, ship body 
paint emissivity, sun elevation angle and sky condition on 
the ship body temperature. 

For both horizontal and vertical faces higher paint 
emissivity and ship internal temperature result in higher 
ship body temperature. Higher internal temperature corre- 
sponds to less ship body heat losses towards the internal 
compartments. Moreover higher paint emissivity (absorp- 
tivity for gray paint) corresponds in higher solar energy 
absorption and, for the 300K ship body temperature level, 
the dominant mechanism of heat transfer is that of convec- 
tion. Hence the increase of the paint emissivity does not 
correspond to a significant increase of the energy radiated 
by the ship body, and the net result is an increase of the 
ship body temperature. 

The sky condition appears to be significant as well. The 
more clear it is the more, solar energy there is incident on 
the ship, resulting in an increase of its body temperature. 

As far as the atmospheric profile is concerned, the 
Tropical Atmosphere results in higher body temperature than 
that of Midlatitude Summer since the ambient temperature is 
higher and the ship body heat losses towards the environment 
are less. 

The sun elevation angle above the horizon affects 
vertical and horizontal faces in different ways. We know 
that the incident solar energy at sea level is an increasing 
function of the elevation angle. For an horizontal face the 
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temperature appears to be an increasing function of the 
elevation angle, the sine of which is used for the calcula- 
tion of the incident angle effect. For vertical faces the 
cosine of the elevation angle is used and the temperature 
appears a maximum in the vicinity of 50° . The explanation is 
that, up to this value, the increasing effect of the solar 
energy dominates the decreasing one of the cosine function. 
After the value of about 50° the first derivative of the 
solar energy function becomes small and the cosine behavior 
dominates so that the net result is a decreasing temperature 
function from 50° to 90°. 

As we see in the ship tracked area geometric analysis is 
constituted from both vertical and horizontal faces, but the 
vertical ones have much larger projected surface areas. 
Therefore the radiant energy from the tracked area is due 
mainly to that of the vertical faces and as we will see from 
the S/N ratio results they appear maximum at a value of 
elevation angle close to 50°. 

B. DETECTOR S/N RATIO 

Figures 5.5 to 5.22 show the effects of the sky condi- 
tion, sun elevation angle and missile flight altitude on the 
detector S/N ratio. 

In general we observe that the variation of the S/N 
ratio is similar to that of the ship body temperature. That 
means that the radiant energy from the ship body due to its 
temperature is the main thermal source for the detector. At 
this point, it is useful to discuss the influence of the 
stack exit plane thermal radiation. From our calculation we 
know that the radiant intensity of the CO in the stack exit 
plane has a magnitude of the order of 10 fc (W/cm 2 Sr) and that 
of the ship body of the order of 10 3 to 10 2 . Moreover the 
stack exit plane area is much smaller than the total tracked 
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area and the stack exit radiant flux is significantly 
smaller than this of the ship body tracked surface. 

Sky condition and the paint emissivity affect S/N ratio 
similarly to the ship body temperature. The more clear is 
the sky and the higher is the paint emissivity the higher is 
the S/N ratio. 

The atmospheric profile affects S/N in a different way. 
While the Tropical Atmosphere gives higher ship body temper- 
ature, it results in lower S/N ratios comparing with the 
Midlatitude Summer. The explanation for that is that the 
increase of the ship body radiance, due to its temperature 
increase, is less than the background noise increase (higher 
sea temperature). Moreover the atmospheric transmittance in 
the 5um window for the Tropical Atmosphere is about 10 % less 
than this for Midlatitude Summer. 

The missile flight altitude appears to be important as a 
parameter. The higher it is the lower is the S/N ratio and 
this decrease is higher between 50 to 70 (m) than between 70 
to 100 (m) . The explanation is that although higher missile 
flight altitude corresponds to lower background area (sea), 
the path between missile head and ship becomes larger and 
its square appears in the denominator of the signal voltage 
formula. Therefore the net result is to give lower values of 
S/N ratio. 

The influence of the sun elevation angle above the 
horizon is similar to that on the ship body temperature and 
it has been discussed earlier. 

C. RECOMMENDATIONS 

Based on the relative influence of each parameter which 
has been varied in our analysis we recommend: 

From the ship design point of view, lower internal 
temperature or, when it is not feasible, better wall 
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insulation, especially in the area where significant heat 
generation exist. In addition to that, paint of low emis- 
sivity and absortivity is desirable. 



From the missile point of view we recommend lower flight 
altitude. At this point we must notice that, lowering the 
flight altitude, we will reach a point where the background 
area (sea) will be very large, resulting in higher back- 
ground noise and the coast or even the sun itself could 
constitute background noise for the detector. Therefore the 
flight altitude must be determined, in a certain design, 
accounting for all these conflicting parameters. Another 
solution could be the use of higher flight altitude, that 
assures us low background area, and to improve the detector 
optics design that will balance the loss of the signal 
voltage due to the greater path length. 
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APPENDIX A 



TABLES FOR SHIP BODY TEMPERATURE 
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Side of Instantaneous Field of View in Meters 
for Flight Altitude= 50 m 
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TABLE 14 

Path in Meters for Flight Altitude= 50 
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TABLE 16 

Path in Meters for Flight Altitude= 60 
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Path in Meters for Flight Altitude-70 



o 


CMCOLOPOtOatiOPOPOioatiOPOMftnOtOPOtd'vOOO 


o 


po po to vo on r-i po to o po mo o> po to o oo csj 


PHI = 5 . 


cm ct to po O <r cm ct to po >— i co to cm O to cm ct 
to to to to to to LO <r <r Ml - MT PO PO PO PO CM CM CM «— 1 1 — i 


.50 


to ct po o cm to «-h r-i cm to to oo po i—i po to r—i o at 
cm to o to o to o to cm oo <r «—i oo to po rH at to to po 


<r 

II 

H 

6 X 

PL. 

3 


<T O po o to po at to cm at to cm at to po at to po o 
io to <t vj - po po cm cm cm «— i r-i i— i o O O at at at at co 
co oo co oo oo oo co oo co co co co oo co co co 


II O 

<u o 


to to •— i cm <r oo <j" <r oo to po o at co i— i r>» to 

r^cMr^PoattocMatr^toPOCMi— 1 «— ii— ir-Hi— icM<J"tnr^ 


Aitituc 
PHI =4 . 


m i — 1 to cm po at <r o to CMOO <r O to CM oo^r O to CM 
to tn<r <r po po cm cm cm i— 1 1— i O O O at at co oo oo 
at at at at at at at at at at at at at at co oo oo oo oo oo oo 


u o 

n to 

hn * 


at to «— i O at tn o at cm<T oo at to to po m o i— i at oo 

to co «— i -d - r— i to cm «— i oo to m tn to at 


uu 

-H PO 

1 — 1 II 

Pn M 

K 

M p-i 

o 


<roopor^r— itootnat<ratpocopooopooopooomco 
oo to to tn m <r po po cm cm i— i «— i O O at at oo oo 

OOOOOOOOOOOOOOOOatatatatat 


rs f 
.00 


tn o at at i— i »— i po to -cr O cm to at uo to cm 

at i — i un oo po o to <r on cm cm po m oo t—i to O vn cm at r>. 


n Metei 
PHI = 3 . 


pn to oo O po to oo i—i-cr o m to at pn to o pn o <r 
lomJ - on pn cm >— i O O at oo oo hx to un uomT mT pn cm cm r— i 


:h i 
.50 


cm pn r-. pn i— i po <r to o pn cn at o oo i— i oo cm o cm at 
at o cm to cm at oo i— i un o to m po <r to oo cm cm 


Pat 
PHI = 2 . 


untnMrpnpncMCMCMCMPnpn<r<fun to r^-oo at i— icm<T 

oo to tn<t po cm i — i O at oo to tn po cm«— i »— i o at 
<r<r<t<t<r<r<tMTMtPOPOPOPOpnpopopnpnpnpocM 


o 

o 


O oo at i— i<r o cm to un i— i to cm to oo oo pn oo 

at oo o cm i— i po oo un un to i-H un to at pn o at at 


PHI =2 . 


pn to o to at <r at <r at tn i— i <r o <r r-i at cm 

cm o at lpkt cm i—i at oo un<j- pn t—i o at to m <r 

cocor^r^r^r^r^r^tototototototototnununinun 


C 

H 

-J 


O CMMT to 00 O CM<J* to oo o cmmt to oo o cmmt to 00 o 
cm cm cm cm cm pn m po po cn <r <r <r <r m un to to 


W 

Q 


ooooooooooooooooooooo 




ooooooooooooooooooooo 



109 



O oin<f cMor^MTOvDcMr^cMvooMrr^ocMMtvDr^ 
O »— ir^coaMn Ovo cm rococo MTOMnomP ivo»— I vd »— 1 



in MDMDr^i^coc^cTNOOj-^»-^cMcMco<r'4"ininvovor^ 



X 

P 



O cNr^^-iunooi— imLnr^r^cor^f^vrxtcNONvDma>tn 

lo oo<ti— ir^ooov^c^Joo^OvOCNjoOvtoio*— ir^cNjoo 



w 


Mi- 


U 


ll 


<u 


H 


u 


X 


<v 


P 






a 


O 


•H 


o 


.5 6 


mt 


<v 


ii 


•HO 


H 


>co 


X 




p 


M-i || 




O <U 




-a 


o 


XJ D 


in 


i— 1 4J 


• 


aj-H 


CO 


H4J 


ii 


1 


H 




X 


w 


p 


a 4-J 




ox 




0 ) QjO 


o 


C'H 


o 


cdi— i 


• 


4->P 


CO 


c 


II 


cd H 


H 


4-J O 


X 


W4H 


p 


C 




M 


o 


m 


m 


O 


CM 


<u 


II 


X3 


H 


•H 


X 


CO 


P 




o 




o 




CM 




II 




H 




X 




p 



vDNoocoa>oOr-iHcNjmnNd"Nd"iovDM)fNNcx)co 



co to vo co vo m cm on vd cm r*^ cm vd o co m oo on 
vd co o <r i— i co in c\j co to cm oo m i— i i oo<t o o cm oo 



oo o' on o »— i »— i cm co co <T in m vo co ctn <tn o O 



io vo<f i — i co m oo cnj vo oo O i — i cnj i — i O a> vd ro 

vO<r cm O oo o <r •— i a> vo vT t— i oo vo m O r^>4* O r^<r 



oo crs o i— 1 1— i cm co <T <? m vo oo o O O »— i cm cm co 



0>4* vo in co o ino covocooo oor^in cm comToo 

o o oo vo m^r co j— i o oo vo < f cm o oo vd < r ■— i a> vo 



O o r-^ cnj m-<r lo vo oo oo o r— i cm CNjm<r lo to vo 

■Ii— 1 1— 1 1— I r-l CM CM CM CM CM CM CM CM CM 



vO LO CM 00 CM tO vO vO lO CM 00 CM to 00 00 00 vO CO CJN 00 

oo a> o 0 1 ~ i — 1 1 — 1 1 — i r— 1 1— i o O a> co vd in-vl" cm r — i o> 



i— i cm^ lh vo oo o i — i cm m cn<r m vo oo a> o O 

i“l i— *1 r— 1 1— I r— H rH r— I t— - 1 CM CM CM CM CM CM CM CM CM CM CM CO CO 



in vo in i— i in in i— i in vo co oo «— i i cm cm r— i oo co 
in oo >— i mt co o cm <r m vo oo on a> o o o o a> a> 



m oo o cm co >4" ^n vo oo a> o cm cn<r «n in vo 

T— I T— I r— 1 1— I T— I CM CM CM CM CM CM CM CM CM CO CO CO CO CO CO CO 



<J OCMMf vOOOOCMMf vOOOOCMmT vOOOOCMMf vOOOO 
H cMcMcMcvjcMcococococo<r<T<f MTMfinininininvo 
hJ 

W ooooooooooooooooooooo 
Q 



ooooooooooooooooooooo 



TABLE 20 

Path in Meters for Flight Altitude= 80 
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Path in Meters for Flight Altitude=90 
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Side of Instantaneous Field of View in Meters 
for Flight Altitude=100 in 
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TABLE 24 

Path in Meters for Flight Altitude=100 in 
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APPENDIX C 



TABLES WITH DATA USED FOR LOWTRAN 6 
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TABLE 25 

Data Used for Month October 



CARD 1 



MODEL 


2 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


2 : 


M2 


2 : 


M3 


2 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


2 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IVSA 


0 : 


VIS 


10 : 


WSS 


2.09: 


WHH 


2.09: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


; 


RANGE 


From 


BETA 


— : 


RO 


: 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 : 



Midlatitude Summer 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Midlatitude Summer T and P Profile 

Midlatitude Summer Water Vapor Profile 

Midlatitude Summer Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Fall - Winter 

No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range(Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/s) 
Rain Rate (mm/K) 



Section G of Chapter II : Initial Altitude 
Table 22 : Final Altitude 

Not Necessary to Define 
Table 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



: Initial Frequency 
: Final Frequency 
Frequency Increment 
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TABLE 27 

Data Used for Month November 



CARD 1 



MODEL 


2 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


2 : 


M2 


2 : 


M3 


2 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


2 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IVSA 


0 : 


VIS 


10 : 


WSS 


2.09: 


WHH 


2.09: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


Form 


ANGLE 


— : 


RANGE 


From 


BETA 


— ; 


RO 




LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Midlatitude Summer 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Midlatitude Summer T and P Profile 

Midlatitude Summer Water Vapor Profile 

Midlatitude Summer Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Fall - Winter 

No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 
Current Wind Speed(m/s) 

24 hours Average Wind S 
Rain Rate (mm/n) 



peed (m/ s ) 



section G of Chapetr II : Initial Altitude 
Cable 22 : Final Altitude 

Not' Necessary to Define 
Cable 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 29 

Data Used for Month December 



CARD 1 



MODEL 


2 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


2 : 


M2 


2 : 


M3 


2 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


2 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IV S A 


0 : 


VIS 


10 : 


WSS 


2.05: 


WHH 


2.05: 


RAINRT 


0 : 
CARD 3 


HI 


F rom 


H2 


From 


ANGLE 


— : 


RANGE 


F rom 


BETA 


— : 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Midlatitude Summer 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Midlatitude Summer T and P Profile 

Midl'atitude Summer Water Vapor Profile 

Midlatitude Summer Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Fall - Winter 

No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range 
Current Wind Speed 
24 hours Average Wind Speed (m/s) 
Rain Rate (mm/n) 



section G of Chapter II : Initial Altitude 
"able 22 : Final Altitude 

Not Necessary to Define 
[■able 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 



Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 30 

Average Transmittances for Month December 
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Total Average of the Month = 0.4704 



TABLE 31 

Data Used for Month January 



CARD 1 



MODEL 


2 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


2 : 


M2 


2 : 


M3 


2 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


2 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IVSA 


0 : 


VIS 


15 : 


WSS 


2.06: 


WHH 


2.06: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


— : 


RANGE 


F rom 


BETA 


— : 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Midlatitude Summer 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Midlatitude Summer T and P Profile 

Midlatitude Summer Water Vapor Profile 

Midlatitude Summer Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Fall - Winter 

No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range 
Current Wind Speed 
24 hours Average Wind Speed (m/s) 
Rain Rate (mm/n) 



Section G of Chapter II : Initial Altitude 
"able 22 : Final Altitude 

Not Necessary to Define 
rable 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 33 

Data Used for Month February 



CARD 1 



MODEL 


2 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


2 : 


M2 


2 : 


M3 


2 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


2 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IVSA 


0 : 


VIS 


15 : 


WSS 


2.22: 


WHH 


2.22: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


; 


RANGE 


From 


BETA 


— ; 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Midlatiitude Summer 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Midlatitude Summer T and P Profile 

Midlatitude Summer Water Vapor Profile 

Midlatitude Summer Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Fall - Winter 

No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/s) 
Rain Rate (mm/h) 



section G of Chapter II : Initial Altitude 
[■able 22 : Final Altitude 

Not' Necessary to Define 
’able 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 34 

Average Transmittances for Month February 
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Total Average of the Month = 0.4710 



TABLE 35 

Data Used for Month March 



CARD 1 



MODEL 


2 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


2 : 


M2 


2 : 


M3 


2 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


2 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IVSA 


0 : 


VIS 


15 : 


wss 


2.13: 


WHH 


2.13: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


— ; 


RANGE 


From 


BETA 


— : 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Midlatitude Summer 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Midlatitude Summer T and P Profile 

Midlatitude Summer Water Vapor Profile 

Midlatitude Summer Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Fall - Winter 

No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/ s ) 
Rain Rate (mm/n) 



Section G of Chapter II : Initial Altitude 
Table 22 : Final Altitude 

Not Necessary to Define 
Table 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 36 

Average Transmittances for Month March 
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Total Average of the Month 



TABLE 37 

Data Used for Month April 



CARD 1 



MODEL 


2 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


2 : 


M2 


2 : 


M3 


2 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


2 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IVSA 


0 : 


VIS 


15 : 


WSS 


2.12: 


WHH 


2.12: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


: 


RANGE 


From 


BETA 


— ; 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Midlatitude Summer 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Midlatitude Summer T and P Profile 

Midlatitude Summer Water Vapor Profile 

Midlatitude Summer Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Fall - Winter 

No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 
Current Wind Speed (m/s) 

24 hours Average Wind Speed 
Rain Rate (mm/n) 



(m/s ) 



?able 22 : Final Altitude 

Not Necessary to Define 
rable 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 



Initial Altitude 
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TABLE 39 

Data Used for Month May 



CARD 1 



MODEL 


1 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


1 : 


M2 


1 : 


M3 


1 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


1 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IV S A 


0 : 


VIS 


15 : 


WSS 


2.02: 


WHH 


2.02: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


— : 


RANGE 


From 


BETA 


— ; 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Tropical Atmosphere 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Tropical Temperature and Presure Profile 

Tropical Water Vapor Profile 

Tropical Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Spring - Summer 
No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/s) 
Rain Rate (mm/n) 



Section G of Chapter II : Initial Altitude 
^able 22 : Final Altitude 

Not Necessary to Define 
’able 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 41 

Data Used for Month June 



CARD 1 



MODEL 


1 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


1 : 


M2 


1 : 


M3 


1 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


1 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IVSA 


0 : 


VIS 


12 : 


WSS 


2.15: 


WHH 


2.15: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


— : 


RANGE 


From 


BETA 


— : 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Tropical Atmosphere 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Tropical Temperature and Presure Profile 

Tropical Water Vapor Profile 

Tropical Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Spring - Summer 
No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/s) 
Rain Rate (mm/n) 



section G of Chapter II : Initial Altitude 
iable 22 : Final Altitude 

Not' Necessary to Define 
[■able 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 42 

Average Transmittances for Month June 
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Total Average of the Month = 0.4262 



TABLE 43 

Data Used for Month July 



CARD 1 



MODEL 


1 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


1 : 


M2 


1 : 


M3 


1 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


1 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IV S A 


0 : 


VIS 


12 : 


WSS 


2.15: 


WHH 


2.15: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


— : 


RANGE 


F rom 


BETA 


; 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Tropical Atmosphere 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Tropical Temperature and Presure Profile 

Tropical Water Vapor Profile 

Tropical Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Spring - Summer 
No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/s) 
Rain Rate (mm/n) 



Section G of Chapter II : Initial Altitude 
?able 22 : Final Altitude 

Not Necessary to Define 
?able 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 



Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 45 

Data Used for Month August 



CARD 1 



MODEL 


1 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


1 : 


M2 


1 : 


M3 


1 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


1 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IVSA 


0 : 


VIS 


10 : 


WSS 


2.20: 


WHH 


2.20: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


— ; 


RANGE 


From 


BETA 


— ; 


RO 




LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Tropical Atmosphere 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Tropical Temperature and Presure Profile 

Tropical Water Vapor Profile 

Tropical Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Spring - Summer 
No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/s) 
Rain Rate (mm/n) 



Section G of Chapter II : Initial Altitude 
Sable 22 : Final Altitude 

Not Necessary to Define 
Sable 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 



137 




138 



TABLE 47 

Data Used for Month September 



CARD 1 



MODEL 


1 : 


ITYPE 


2 : 


IEMSCT 


0 : 


Ml 


1 : 


M2 


1 : 


M3 


1 : 


IM 


0 : 


NOPRT 


0 : 


TBOUND 


0 : 


SALB 


0 : 
CARD 2 


IHAZE 


3 : 


ISEASN 


1 : 


IVULCN 


0 : 


ICSTL 


5 : 


ICIR 


0 : 


IV S A 


0 : 


VIS 


8 : 


WSS 


2.17: 


WHH 


2.17: 


RAINRT 


0 : 
CARD 3 


HI 


From 


H2 


From 


ANGLE 


— ; 


RANGE 


From 


BETA 


— 


RO 


— 


LEN 


0 : 
CARD 4 


VI 


2255 


V2 


2000 


DV 


10 



Tropical Atmosphere 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 

Tropical Temperature and Presure Profile 

Tropical Water Vapor Profile 

Tropical Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 



Navy Maritime Extinction 
Spring - Summer 
No Stratospheric Background 
Medium Continental Influence 
No Cirrus 
Not Used 

Meteorological Range (Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/s) 
Rain Rate (mm/n) 



iection G of Chapter II : Initial Altitude 
[■able 22 : Final Altitude 

Not Necessary to Define 
Cable 22 : Path Length 

Not Necessary to Define 
Not Necessary to Define 
Normal Operation 



Initial Frequency 
Final Frequency 
Frequency Increment 
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APPENDIX D 



PLOTS OF S/N RATIO 
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UNOBSCURED SUN 
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Figure D.2 S/N Ratio for Midlatitude Summer, £*=0.80 
and Flight Altitude 100m 
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APPENDIX E 



COMPUTER PROGRAMS 
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on n nnnnnn 



$ JO B 



PROGRAM SHIP1 

» * * 



CALCULATION OF SICE UF INSTANTANEOUS FIELD OF VIEW 

DIMENSION ti (b) ,PHI(7) , DELTA (21) ,SIDE(o,7,21) , PATH ( t> # 7 r 2 1 ) 

DATA (H ( I) ,I=1,6)/50.,60.,7C.,80.,90.,100./ 

DATA (PHI (K) ,*=1,7)/2.0,2.5< 3.0,3.5 c 4.C.4.‘j f 5.0/ 

DATA (DEL: A (L) ,L = 1 , 21)/. w,. *2, - <<4,. zb, . zd, . j, . 32, • 34, . Jb, .38, 
1 • 4 # .42 # .44,.4b,.4<J,.5,.52,.54 # .5c,.58,.o/ 

DO 50 1=1,6 
DO 45 K= 1 , 7 
DO 40 L= 1 , 2 1 

TH= (90.-PHI (K) -DELTA (L))*3.141592/1o0. 

B = U I) * I A N l 1 H) 

C=(H (1) *H (I) ♦u*b> **.5 
' “ " - - - ■ >92/180- 



C= ( H (1) * H (I) ♦ u*b> **. 
D=DE-TA ll) *3. 1'4 159 2/ 
SIDE (I, X, Lj =2. *C*SI t 



40 

45 

39 



38 



34 

35 



, *C*SIN (D) 

PATH ( I , K , L) =SIDL (i,K,L)/ (2. *T A N (D) ) 

CONTINUE 
CONTINUE 
WRITE(o f J9) ri ( I) 

FUEulAT (* 1 ' 4< ■ SICE OF I MST A NT* N EGOS FIELD OF VIEW IN 3 ETERS 1 // 
1bX,'FOn FLIGHT H £ I GHT= ' F 1 • 0 , 1 X , • 3 El E 22 ' // ) 

WRITE (6, 3d) (PHI (XI , k- 1 f 7) 

FORMAT (2 X , ' DELTA' , 3 X, 7 ( ' PH I = 1 ,14.2,21)/) 

CO Jb L = 1 #21 

WHITE (b. 34) DELTA (L) , (81CE(I,K,L) ,*=1,7) 

^ IQ* ,2X,F4.j/U, 7f 10.2) 



FORilAt ( 
CONTI NUi 



W PI IE lb . 49) H ( I) 

49 FORMAT (M • 4X, ' PATH IN JETERS •// 3 X FOR FLIGHT HE IG H X = • , F 4 - J , 

1 iy JlFTFR ^ ■ y / \ 



48 



46 

47 
50 



IX * METERS'//) 

7 (*PH 

WRITE (b f 4b) DELIA (L) , ( t- A T H ( I , K , L ) , N = I , 7 ) 
FORNAf j 5 0', 2X, ?4. 2, )X, 7F U.2) 



WHITE { o , 4 3 ) l?HI(K) , K = 1 . 7 ) 

FORriAi Ul, • DELTA 1 , JX , 7 ( U'HI = » , ?4. 2 , 2 X ) /) 
CO 47 L = 1 ,21 



IENT8X 



CCIUIHU& 
CONTINUE 
ST C P 
END 
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luU'J UU U U UUU'JU uuuuu uuu uu 



SJOB 



PROG ii A M 



jIII? THERMAL 31GDAIJRE 

l *********** * 



DIMENSION T (2 , 2, 3,91) , £I& (2) ,13(2-91) -ESN (2, 3,91) , C 
DIMENSION HDL (2,2, 3,31) ,HDLR (2, 3.91) . fiDLT i 2 , 2 , 3 , 9 1) 
DIMENSION LiI13,9l),Rl^(J,91),jl3(3,3l),aI4(J,91) 



0 (91) 



DIMENSION U 1 3 ( 3 , 3 1 ) ,317(3,91) 
DIMENSION VS (J ,9 1) , JNh ( 3, J 1) 



DATA (TIN ( J ) ,J = 1,2) / 29 3. 15,303. 15/ 
DATA EIN f ?I/. 44, j. 1*1552/ 

DATA TO, I HANS/ 294. 2, J. 471 1/ 

DATA TO, TRAN S/299- 7,0. 4259/ 

DATA EO, AO/0.30,0. 30/ 

DATA EO, AO/O. 35- J. 35/ 




DATA HEIGHT, SIDE, 2 ATS/ UdUU. ,2^16.0, 11 j 372.0/ 
DATA PATHS, AREAS/ IdOobd. , 1 -U^0 30./ 

DATA PAT HS , >* 8 E A S/ 1b4o07., 1316u00./ 

DATA PATHS, An E AS/ 154747. , U 77 6 00./ 

DATA PATHS, AF EAS/ 146194. , UI2300./ 

DATA PATHS, AREAS/ 143 35 0., 1 14 4 200./ 

DATA PATHS , A H E AS/ 1 3 9 b 0 7 . , 1 llbOJO./ 

DATA HCONV/5o. 750/ 

DATA APERT, ThGPI, RESDON/1 1j. 1,.J1.2£*04/ 

DATA tuCL, XTOT.EN, L’l MZr , 3 ET t C. 1/ 4 . 24 - 500C 3. , 1 . 
DATA Cl, CZ,S6/3. 74156-12, 1 . 4 J H d , b.o/E-08/ 



,. 01 , 56 * 10 / 



DATA ( O (L) , L=1,91,5)/0.,D.,ig.,15.,20.,25.,30.,J5.,40.,45.,50., 
55. , u d . , b 5 . ,70. ,7b. ,o0. < u 5 . , 90./ 

DATA (£S (1,1). , L= 1 , 3 1 ,5)/b. / ,4d.b,9D.5, 147.0, 19 5. 7, 25b. 7, 319. 9, 

1 39 1.3, 47 -9. 9, oo4. J , b J o . 7 , 7 0 3 . 1 , 7 o 3 . 7 , d 1 j . o , 3 5 4 . C, 379.9,093. 1 , 

1 893. 1,930.0/ 

DATA (ES (2.L) , L= i, 91, 5)/2. 7, 1 9 . 7 , j 9 . 1,64.4.9 1.7, 117.4, 140.5, 

1 181. 7,2 19. 1, 25 4. d, 3 0 0. J, j4z. 2 , j 79. 3, 413. 5,444. 1,4b 9.6,493. 4, 

1 5 1 3 . 2 , d 3 2 . 2 / 

DATA (ES( 3.1.) ,L=T .'J 1 , 5) /3 . 5 , 4 . 3 , 9 . 0 , 1 3 . 7 . 1 9 . o t 25. 4 , j2. 3-43.0, 

1 53. 2, ol. 0,o7. 7, 74 . j, 79. 9, 85.1, 39.3,94. 2, 98. 5, 102.1,1 05. 7/ 



C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



THE HE LOW DaTA ARE IN CM OR CM^ 

DATA - 1 , H 1 ,82, H3/ 162 0. .7 23. , 2 10 . , 900./ 

DATA V4,il4,HS, *5, il 7/90 0. , Jcu. , 1SJ. ,o 18. , 12 60 
DATA A4 , a <IoN/->. 0 73E + U5,oo. 47b/ 

DATA HEIGHT, SIDE, PATH/o 00 0.0, 2 1b7. , 1 1^339./ 
DATA PATHS, m RE A 5/ 1b4 307. , 1 jIcc'JO./ 

DATA £U , 1 J/0. 9 J , ^9 15/ 

DATA TRANS/3. 42o9/ 

DATA (WE (N) , 11 = 1 , 43) / 5. 003, 4. Odd, 4.9 7 5,4.963, 
14. 914, 4.302,4.350, 4.678, *4 • d 6 o , 4.354,4.343,4. 
14. 735, 4. 773,4.76^,4. 751, 4. 73 9, 4. 728, 4. 717, 4. 
14. do.*, 4 . ojI, 4. 64 0, 4. b30, 4. ol9,4.o0d, 4. 598,4. 

DATA (HL( M),N = u1,a2)/4. 34b, 4. S35, 4.325,4.315 
14.474,4.404,4.454,4. 444.4. 4 j5/ 

DATA (TRL(N) , 3=1, jJ)/ 0.1 o 77,0. 1oo3, 3. 2163,0. 
10. 2917,3. 3334,0 .36 19, J. 3 79 7,0. ibcy, 0.3902, 3. 
10. 4747, 0. 5732, 3. b2o 0,3. 0427, 0.o4Jd, 0.0^24, 3. 
10.b8do,j.oo17,3.obbD,0.o^o7,0.ojob,3.bbC2,3. 

DATA (TRI(N) - ii = 3 4 5 2) / 0 . 7 3 39 . 3 . 7 2 b □ , 3 . 7 1 3 1 , 3 
10 . d342,J. 4940, J.4od 3, 0.4335,0. 4174, J.350o,J. 
10.2099 ‘ 



0.20 99,3. 14o 6,3. 1 1 4b, 0. 39 3 8/ 

DATA {TSL (N) -3=1, 3 3)/ 0.399 0,3. 112 2, 3. 14 57, 3. 
0.216b, 3.2550, 3.2813, 3.3304,3.3034,3.314 j,0. . 
0.3995,3 . jO 17, 3.552o, 3.5733,3. j74d, 3. 552 7,0.: 
10. b057, 3. 5 97 4, 3. o9zJ,0. 533 5,0.5 74 7,3.5394.3.1 
DATA (TSL ( N ) ,N=J4,52)/0.bd9b,J.bob7,0.o75b,0, 
10. 5113,3.4 73 5, 0.4 460, 3. 4 07 5, 0.4 30 0,3. 3760,0.. 



./ 



4. 350 , 4.938,4.92o, 
031,4.819, 4. 808, 4. 79b, 
706 ,4 . o9b , 4.684,4.673, 
587,4.577,4. 56o, 4. 5b6/ 
,4.305,4. 494, 4.484, 

2070.0. 2070. 3. 2448, 

3852.0. j 794, 0.4242, 
oj43,3.bb25, 3.6331, 
c643, 3.b8o8,3. 7053/ 

. oc99-0.o277, 3.5907, 
3589, 5. J^21,d.2b20, 



1 4 3 o , 0 . 142 1,0. 1727, 

. 3126,0.3035,3.3511, 

10.3995,3.j017,3.552o, 3.5733, 0. j74d,J.5527,0.Sc57, 0.5972,0.6217, 

4 " "" " " " ‘ 1 " <r * ,r ‘ " " " ' ' ' 60 j 3 . 3. b 29^ - 0 . b5u2/ 

.636^,0. 5374.0. 5b41 , 
3469, 3. j 123, 3. 2555, 
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nnn no 



1 

2 

J 



10,2046, J. 1 4 J d , J . 1127,U.U J2 *4/ 



DC 3 K = 1 , ^ 

DC 2 11= 1 , Id 
N N= 5 ♦ 1 1 ♦* 1 
N = 5 ♦ ( 1 1 - 1 ) o 

do i l = *,;:n 

ES ( K , L ) = 23 ( K , N J ) - ( 2 S (a, 
D (L) = D (UN) - (0 (NN)-O (J) ) 
CONTINUE 



:i N) -13 
♦FLOAT 



(X, N ) ) ♦FLOAT ( N N - L ) / 5 . 
( N N - L ) / 2 . 



CONTI NJE 
CONTINUE 

CALCUL AT-iON 0? N C n .1 A i- CC.1FCNENT3 OF ES 



C 

C 

C 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



7 

3 

9 

10 
1 1 

12 

1 J 
14 



15 

16 



DO 6 1=1.2 
DC 5 K= I , J 
DC 4 L=1,*1 
1 H = 0 (L) ♦FI/160. 

Ir (I. -3 - 1) A = JIH < III) 

IF 1. 20. 2 A = C 3 : » ( I’a) 

ESN (I, K, L =CS (K . L) ♦ \ 

BDLli(I,2,L) = (9.J4b-07; ♦ESN ( I , K , L) 

CCNIINUE 

CONTINUE 

CONTINUE 

WHITS (b .7) 

FGRMATI M * ,dX, ' ScHK CCMSIA -A i -3 IN A/*2'///) 

WHITE (o,o) 

FORMAT (IX, *ALi. ANGLE' , UX, 'U 40 jSC’JFED SUN*, 2a, ‘SUN «IIU LIGiil CLCU 
1 DS * , 2 X , 1 3U Wild HEAVY STCi.rt CLOUDS'/) 

DO 10 L = 1 , i 1 

WRITE lb j) 3 (L) , (23 iK,*.) ,K=1 .J) 

FORM AT ( 1 X , F o . 3, oa , iM J. 2 , oX ,f 1j. 2, 14 X, F 10.*) 

CONTINUE 



WRITE (o. 11) 

FOhMAT (* 1 ' # » bCdJAL COMPONENT OF F.XDI.-iNCE DUS TO LEFT EC TED SCLAH EN 
1£KGY*//1X, ' IN ^/(ClI^SF) I k > Till WINDOW C F 5 :i IC HC M E I E I •///) 

W R I T E ( o . 1 c ) 

FOadAi (‘O' Ji, 'IJNOuSCUn^D jHN' ,32, '30; 

1 *3LiN .i IT H HEAVY 3 I C i : *1 CIOUDS', 



« IT 11 L I w»d £ LL C J j ‘ , 2 X , 



2000 



WRIIE (6. U) 

FOKJA1 M X, ' ALTITUDE* , 1 X, J ( ' HCE .-OUIAL* ,2X, ‘ VEnllCAL ‘ , 2X) ) 

WHITE (o,14) 

POEM A I ( 1 a , • ANGLE ‘ , 4 k ,3 ( ‘ t Acr. ‘ , UX, • 5 ACE * ,fcX) //) 

DU 16 L = 1 , 9 1 

WHITE (o, 15) d (L) , (HOLD (1, 1,1) ,1= l,.) , UOL?. (i, 2,1) ,1=1 ,1) , 

I (;]2a.i\(I,J,I),j.= 1,^) 

FORMAT (1X,F4.J,j\,j(2 1u. 4, 1 a , L 1 J . 4, IX)) 

CONTINUE 

CALCULATION UF Still BODY EEMHEHAIUHS 

DC 20 1=1,2 
DC 19 J=1,2 
DO Id :<=i,3 
DO 17 L — I , J 1 
io4=to»*4 . 

TIN4=TIN (.1) * ♦ 4 . 

ESN ( I , K , L) = (../J.) ♦ESN { I , K , LI 

A= ESN ( I , K , L) *AC* , 20*3j*TC4*’EiN*Sfc* , :iN4 

d= ( £0 ♦ c 1 N ) * 3 D 

X K= (A/D) **9. 2 J 

T4=XK**4. 

FI — EG* SB* (T 4 -TC 4 ) ♦ El w *5 u* (14- TIN 4) ♦HCONV ♦ ( X X - 1 0 ) - E 3 N (I, K, L) ♦AC 

X A = X K * F 1 

X4=XA**4. 

F2 = EO*S3 + (X4-T04) fEi.N*Sfc* (X 4 -I ID 4) ♦HCCNV* (XA-TO) -£3N (I,X, L) »AC 



X K 1 = X K - r 1 ♦ c 1/ I 
D I ? = A u S (Xfv1-<K 
X K= X K 1 



F2-F1) 

) 
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nnnr.nn noon noon onnn nor. onon none noon 



17 

10 



If {DIF.GT.000U01) GO 10 2000 
T(I,J,K V L) = X K 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 



21 

22 

23 

24 



25 

2b 

27 



DO 27 J= 1 . 2 
W3IIE (6.21) TIN (J) 

FORMAT (M ' 3X, ’SHIP BOD 7 Tt-MPEH AT URE IN K FOB T I N = ' , ? 7 . 2 , 1 X , 
1 1 K * /// ) 

• HITE (to. 22) 

FORMAT ( * U 1 , 9 X , 1 T J N 0 3S CU 3 E 0 SUN*, JX, 'SUN WITH LIGHT CLCUD',2X, 
1 'GUN WITd HEAVY JiURrt CLOUDS*/) 

WRITE (o. 2 J) 

FORMAT (IX ' ALTITUDE* , IX, 3 ( MiCR 120 NTAI 1 , 2X, 1 VERTICAL • , 2X) ) 

write (q,A) 

FORMAT (U. • dN OLE • , 4 X , 3 (• PACE* ,dX, • FACE* , t>X) //) 

DO 26 L=1,J1 

WHITE (0,^5) J(L) , (T (I,J, 1,L) ,1=1,2) , (T (I,J,2,L) ,1=1,2) , 

1 (T (I,J, 3,1) ,1=1 , ^ ) 

FORMAT ( 1 X , f b . 1,3 (?10.2,lX,PlQ.2, IX) ) 

CONTINUE 

CONTINUE 



CALCULATION OF NORMAL COMPONENT OP SPECTRAL RADIANCE IN THE 
WINDOW Of 5 MICROMETER (I-E. 4.435-5. 0 C C) 



23 



29 

30 

31 

32 



DC 32 1=1,2 
DO 31 J = 1 . 2 
DO 30 K= 1 , J 
DO 2* L=1,31 
SI =0.0 
Sz=0. 0 



XI =02/ (T 
“2 = 0 2/ (I 



i,o,k,l; 
i,j,k,l: 
23 M = 1 t 0 



X 2 

DC 

A 1 = FLu A r 
A2=r LCAT 
b = F L C A T ( 
S1«SU b) 
S 2 = S2 ♦ 



♦4. 435E-04) 
^b.OOOE-O*) 



U) *xi 

M) * X z 
A) **4. 



s i=s ueap'(-a i ) * (a i** j. f3.nmub.Mii5.; 

■ XP (- A 2) * AZ- ' ~ 

CONTINUE 



• ♦* 3. *-3. *A2*A2*c. *Al + b. 



/0 

/B 



HDLI {I,J,.X,L) = C 1 * EC * ( T ( I . J , K , L ) * * 4 . ) * (32-S 1) /(?I*(C2**4.)) 
HDL ( I , J , K , L ) = H D L T (i,J,R,L) ♦H.OX.fi (I,K,L) 

CONTINUE 



CONTINUE 

CONTINUE 

CONTINUE 



CALCULATION OF NORMAL COMPONENT OP SEA RADIANCE IN THE 
WINDOW Of 5 MICROMETER (I.E. 4.43S-5.CGC) 

S1=0.0 
S 2 = 0 . 0 

X 1 =C2/ ITW *4 . 4 35 E- 34 
X2=C2/ (TW+5.000E-34 
DO 123 11=1, ti 
A 1= FLOAT (M) *X 1 
A2 = FLUAT (nj * X 2 
B= PLO A I (M) ♦♦a. 

S 1 = 31 4 EXP (- A 1 ) *(Al*n.f 3 .nlnUb.nlfO.)/B 
S2=S2*EX2 -Ai*A2**3.*3.*A2*Az+6.*A2*o./3 
128 CONTINUE 

HDL W = C 1 * E W * (T W * * 4 . ) * ( S2-S 1) / (PI* (C 2**4.)) 



DO 40 J= 1 , 2 
WRITE (b .3 3) 

33 PORMAI ( '1 ' , IX, • NORMAL COMPONENT OF RADIANCE IK w/ (CM2*S3) '//IX, 
I'DUE TU THERMAL a A C IA I* ION * /) 
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34 

35 

36 

37 



38 

39 

40 



WBITEfo, 34) UN (3) 

FORMA! (zl, ' THE WINuGW OF 5 MIC HCM ETE B PCS TI N = ' , F7 . 2 , 1 X , • K '///) 
W B IT E ( 6 . 3 5 ) ' 

FO BM A I ( • 0 1 , 4 X , 'UNOBSCURED JU N 1 , 3X , 1 SOU WITH LIGHT CLCUD',2X, 

1 'SUN .lid HEAVY STORM CLOUDS'/) 

WHITE fo . 3o) 

FORMA! Ml, ' ALTITUDE' , IX, J (• HORIZONTAL 1 ,2X, ' VERTICAL ' , 2X) ) 

WRITE (6.37) 

FOB MAI ( IX, ' ANGLE ' ,4X, 3 ( ' FACE ' , 81 , 1 FACE' , 6 X ) //) 

DO 39 L= 1 , 9 1 
WHITE (o , 3d) 

FOBMAT (II, to 
CONTINUE 
CONTINUE 



. 1 , 4. X , J ( G 1 J • 4 / 1 I , L 1 O • 4 , 



J, 1 ,L) ^1=1,2) , (HDL(I, J,2,L) ,1 = 1,2), 
H) ) 



WHITE (6 . 138) TW, HDLW 

8 FCSM A I ( * 1' , 2 X , 'SEA RADIANCE FOR T W= 1 , F6 . 2 , 2 X , E 1 5 . 4 ) 

CALCULATION OF SHI? SURFACES AND THEIR CORRESPONDING MEAN ELEVATIONS 
FBOfl THE SEA, INSIDE THE INSTANTANEOUS OETSCIOn FIELD CF VIEW IN THE 
B EG IQ b OF THE HIGH EH RADIAN! INTENSITY* 



A 1 = H 1 * tl 1 

All-- (SIDE-W1) ♦HI 
HM 1=81/2. 

A2=Wl*H2 

A22= (SIDE-4 1) *H2 
H M 2 = H 1 



A3=SID£*H3 

Hfl3=Hl+H3/2. 

HM 4 = H 1 + HJ *■ ri 4/ 2 . 

Ab=Wb*H5 

HM5 = ii 1 +U3*H4*d5 //. 

HM 6 = H 1 ♦ HJ + 114 + iib 

A7=SIDE*H7-PI*W4**4/4. 
EM7 = Hl*tiJ 



CALCULATION OF THE RADIANT INTENSITY IN (»/SH) OF THE 
SURFACES (AUA11) , (A2*A22) , AJ,A4,A5,A7 IN THE WINDOW OF 
5 MlCaOtlETEH , IN THE DxpECIICN SHIP-MISSILE 



DO 42 K 
DO 4 1 
Bl IN 
C 1 1 = 
Cl® ( 
B 1 1 
B I2N 
C2 2 = 
BI2 ( 
RI3N 
C33 = 
C3= ( 
BI3 ( 
H 1 4 N 
CU4 = 
C4= ( 
BI4 
B 1 5 N 
C5b = 
C5 = ( 
Rib { 
C 66 = 
B 1 6 = 
HI 7 N 



s 1 / 3 
L = 1 9 1 

= 1131/(2, 2, X, L) *A1*8CL (2, 1 ,K,L) 
(HEIvjHI-hMl) /PA.IH 
1.-C1 1*C1 1) * * 0 • 5 
K,L) =5I1N*C.1 

= bDL(l,2,X,L) * A 2 * H D L ( 1 , 1,K,L) 
(HE! JllT-HM2) /PATH 
a , L ) = 3I2N*C22 
=HDL(2,1,K,L) M3 
(HEIuHT-HM J) /PATH 
1.-c33*C33j ♦♦0.5 
K,L) =uI33*C3 
= tiiJL U , 2, ri , L) * A4 
(HEIGHT- 8 fl 4 ) / P A I H 
1. -C44*C44) *♦0. S 
K , L ) =RI4N*c4 
= H DL U,2, K £ i.) ♦ A 5 
(riEIGHI-hflO) / PA IH 
i. -C55*C55) ♦♦0. 5 
K , L) = B 15 N 5 
( H£IGHT-rii16) /PATH 
aI6N*C6o 
L10L<1, 1 , K , L ) ♦ A 7 



♦A 1 1 



* Az 2 
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onnn nnn non r.nn 



41 

42 



C77 = ( 3E IGriT- H a7) / 
317 (K,L) =HI7N*C77 
CON I IN UE 
CONTINUE 



pa ia 



CALCULATION OF SHIP SIGNAL VO L I AG 2 



49 

50 



DO 50 K= 1 , 3 
DO 49 L = 1 , J 1 

YS ( K , L) = ( R I 1 ( K , L ) ♦ R I 2 ( K ,L) *3 13 (*. L) +3I4(K,L) >315 (K - L) ♦SI6 + 
RI7 (K,L) ) *TRAN5*h£SFCa*0.d 1 * A ? i fi T/ (FATH*FAIH) 

CONTINUE 
CON f I N U I 



CALCULATION OF NOISE VOLTAGE 



OS I NST =SI DE*i IDE/ (PATH *i? AT H) 

AD = 0:JIMST*FCCL*FOCL 

Oa iOT= (SIDE/ PATH) *ALuG ( (LA In* <TOT) / (PATH- XTCT) ) 
DFaU* = UJTOT/ (2. *UiiINSI*EN*TI.12F) 

YMOPi*aESFON* ( ]A0»CFB£:) ** . 3)/L£T£CT 

VB=dQLW* AREAS* I RAWS *RSS?CW *J.d 1*APE£I/ (?ATHS*P AT H5) 
VNTOT-VNCPI^Va 



CALCU LA1XON OF SIGNAL TC NOISE RATIO 

DO 60 K=1,J 
DO 09 1=1,91 

SNP (K,LJ =20. * A LOG 10 (VS (K ,L) /VNTOl) 

59 CONTINUE 

60 CONTINUE 



TAJULATION OF THE DETECTOR SIGNAL TO NOISE RATIO 
yaiiE{6,70) 

70 FO&aAi (' 1 • ,aX, • DETECTOR SIGNAL TO NOISE RATIO* 

1/8jC,*FOR f LI jH I HEIG HT= 50 MEIERS AND SICE=21.Jd M El E ES ' //) 

fa fill E (b r 7 1 ) 

71 FORMAT (»0* , ‘ALT. ANGLE 1 ,31, ‘UufiflSCUFED SUN*,JX,*SUN ^ 1 1 H LIGHT CLOD 
1D*,JX,'S0N WITH HEnVY SIORM C1CUDS'//) 

DO 73 L= 1-01 
WRITE (b, 7 

72 FOB H AT ( U 

73 CONTINUE 
STOP 
END 

SENTBY 



2 ) 



Fd 



D (L\ f (SNP. (K,L) , K= 1 , 3) 

J, 7 a, FI 0.3, 1 3 X , F 1 0 . 3) 
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